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We consider a possible detector-efficiency loophole in experiments that detect entanglement via the 
local measurement of witness operators. Here, only local properties of the detectors are known. We 
derive a general threshold for the detector efficiencies which guarantees that a negative expectation 
value of a witness is due to entanglement, rather than to erroneous detectors. This threshold 
depends on the local decomposition of the witness and its measured expectation value. For two- 
qubit witnesses we find the local operator decomposition that is optimal with respect to closing the 
loophole. 
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I. INTRODUCTION 

Entanglement is the central theme in quantum infor- 
mation processing. It allows to design faster algorithms 
than classically, to communicate in a secure way, or to 
perform protocols that have no classical analogue. En- 
tangled states of few particles (e.g. photons, ions) can be 
routinely created in experiments, and their entanglement 
can be confirmed using state tomography, Bell inequal- 
ities or entanglement witnesses. All of these tools are 
well-established methods for the detection of entangle- 
ment. But can one be sure that they give a confirma- 
tive answer even when realistic, i.e. erroneous detectors 
are used? Here, we will introduce and study a loophole- 
problem for the detection of entanglement via witness 
operators. 

Loophole-problems have been widely discussed in the 
context of ruling out local hidden variable (LHV) mod- 
els, by measuring a violation of certain inequalities, as 
suggested in the seminal work of J.S. Bell in 1964 
Many experiments have been carried out along that line 
Q , but all of them so far suffer from the locality loop- 
hole (i.e. no causal separation of the detectors) and/or 
the detection loophole (i.e. low detector efficiency). As a 
consequence of a loophole, quantum correlations are also 
explainable by LHV theories @, 0] ■ 

In this paper we discuss a possible detection loophole 
for experiments that measure entanglement witnesses. 
Here, the goal is not to prove the completeness of quan- 
tum mechanics (as in Bell experiments), but, assuming 
the correctness of quantum mechanics, to prove the ex- 
istence of entanglement in a given state. One advan- 
tage of witness operators is that they require only few 
local measurements to detect entanglement Q; global 
measurements are experimentally not easily accessible 
at present. A local projection measurement with real- 
istic imperfect detectors (in the computational basis, for 
qubits and isotropic noise) can be described by the fol- 
lowing positive operator valued measurement (POVM): 
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F = £ |0) (0| , F 1 = £ |1) (1| , and F 2 = (1 - £)1, where £ 
is the efficiency of the detector. However, in general the 
global properties of the detectors are not fully charac- 
terised, e.g. there may exist correlations between POVM 
elements of different detectors. Provided that only local 
detector properties are given, what are the conditions for 
being nevertheless able to prove the existence of entan- 
glement without doubt? 



An entanglement witness W is a Hermitian operator 
that fulfils tr (Wp s ) > for all separable n-partite states 

p s = J2iPi^>2 ){ipv | 0], where the index v num- 
bers the subsystem, the probabilities pi are real and non- 
negative, J2iPi — L an d tr (W/9 e ) < for at least one 
entangled state p e @, Q. Throughout this paper, we will 
use without loss of generality normalized witnesses, i.e. 
tr(W) = 1. Our goal is to ensure that a negative mea- 
sured expectation value (W) m < is really due to the 
state being entangled, rather than to imperfect detectors. 
The following line of arguments also holds for specialized 
witnesses which are constructed such that, e.g., they de- 
tect only genuine multi-partite entanglement [9J] or states 
prohibiting LHV models (loj . 



This paper is organized as follows. After introducing 
the local decomposition of entanglement witnesses, we 
study the effect of lost events as well as additional events 
on the experimental expectation value of W. Here, we 
use the worst case approach to derive inequalities which 
need to be fulfilled to ensure entanglement of the given 
state. The parameters in these inequalities are the mea- 
sured expectation value of the witness, the detection effi- 
ciencies, and the coefficients for the local decomposition 
of the witness. We show for the two-qubit case how to 
optimize the local operator decomposition of the witness, 
such that the detection efficiency which is needed to close 
the loophole is minimized. Some recent experiments mea- 
suring witness operators are presented, to demonstrate 
that the detection loophole is a problem in current ex- 
periments. 
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II. DETECTION LOOPHOLE FOR 
ENTANGLEMENT WITNESSES 

Any witness for an n-partite quantum state in d = 
T\™d u dimensions can be decomposed in a local operator 



basis, i.e. an n-fold tensor product of operators 



(i) 



w 



i=0 



(1) 



— (i) 

where the coefficients c, are real. Each operator o v is 
traceless or the identity and corresponds to the ith local 
setting for the party number v. In this expansion, we 
include implicitly also the local identity operators which 
do not need to be measured. The number of terms in eq. 
([1]) depends on the decomposition, i.e. on the choice of 

operators \ A straightforward, but not necessarily op- 
timal choice (concerning the needed detector efficiency) 
are the d\ — \ Hcrmitian generators ct„ of SU(<i„) and the 
respective identity operators. 

In the following, we will use a simpler notation, namely 



(2) 



a=l 



where S a stands for one term from the local expansion 
(fT|). Here, we exclude the identity 1 (acting on the total 
space) from the sum over a, because it does not have to 
be measured and therefore has a special role. 

We will now investigate one local measurement setting 
described by c a S a , and drop the index a for convenience. 
The measured expectation value c {S) is given by 



■(S)m = 



- £-») Ai 



N 



N + e+ - e_ 



(3) 



where Ai is the «th eigenvalue of S, the number of mea- 
sured events for the ith outcome is denoted as 7ij , and N 
is the total number of measured events for that setting. 
In the second part of eq. (J3j) we expressed this expecta- 
tion value as a sum of the ideal number of events, denoted 
as hi (i.e. for perfect detectors), the additional events £+i 
(e.g. dark counts) and the lost events e_i for the ith out- 
come. The total ideal number of events for such a setting 
is denoted as N, and the total number of additional/lost 
events as e±. We have TV = hi and e± = J2i £ ±i- 

The experimental data usually gives no information 
about the number of errors s±i for a specific measurement 
outcome i. Only the detection imperfections are known, 
namely the detection efficiency ("lost events efficiency"): 



N-e- 



V- = 



N 



and the "additional events efficiency" : 

N 



N + e+' 



(4) 
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where < r}± < 1 holds. Here, r)± denotes the global de- 
tection efficiency for a given measurement setting. In the 
following, we assume that rj± is the same for all settings. 
(Other cases can be included by indexing rj± with a.) 

Usually, one makes the fair-sampling assumption about 
the statistical distribution of the unknown errors e±i, i.e. 
one assumes the same statistical distribution for detected 
and lost events; the additional events are assumed to have 
a flat distribution. Here, we will give up this assumption 
and will consider the worst case, where both lost and 
additional events contribute such that the expectation 
value of the witness is shifted towards negative values. 
We point out that in order to reach this worst case sce- 
nario, it is already sufficient that the global POVM ele- 
ments exhibit certain classical correlations (while being 
compatible with the local measurement operators) [llj . 

The worst case is equivalent to finding the lowest pos- 
sible c(S) m , which is achieved by minimizing the con- 
tribution of the additional events, and maximizing the 
contribution of the lost events in eq. This mini- 

mization/maximization can be easily shown to have the 
formmin £+! cJ^e+jAj = e + fl + andmax E _ t cJ2i e -i^i = 
with 

fi± = c (9(±c)A min + 9(Tc)A max ) , (6) 

where 0(x) denotes the Heaviside function, and A r 
is the minimal/maximal eigenvalue of S. 

Inserting this into eq. (J3J), one finds the following worst 
case estimate for the measured expectation value c (S) m : 



mm / max 



where k 



c(S) t -n + (i 



N(N + e+ 



V+' 



— + V- 
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Here, we have introduced the notation c {S) t for the true 
expectation value (without any errors). 

Using (W) m/t = c + J2 a c a ( S <*) m /t and e q- 0. with 
isotropic detection efficiencies, we can express the mea- 
sured witness expectation value {W) m as a function of 
the true one (W) t , for the worst case. To close the loop- 
hole it is necessary to ensure that (W) t < 0. This leads 
to a condition for the maximal {W) m that depends on 
the decomposition of W and the efficiencies r]±: 



(W) r 



< co(l - k) - 



1- — ) + n„_(i-»7_) ),(8) 



where we have re-introduced the summation index a. 

/From now on we want to focus on the case where the 
subsystems are two-dimensional, i.e. qubits. For qubits 
the measurement operators S a are chosen to be tensor 
products of Pauli operators with eigenvalues A = ±1. 
This simplifies f2± to £l a ± = T \c a \, and eq. |(5J) reads for 
qubits 



(W) m < c 
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FIG. 1: The contours correspond to the maximal {W) m of a 
multipartite qubit witness which still ensures entanglement, 
where a decomposition with Co + *^2 a \c a \ = 1 was assumed, 
which is e.g. achieved for an optimal two qubit Bell state 
witness Wu,+ (see text). 



such that we get the Schmidt operator decomposition 

3 

W 3 = ^Sid A i®a B i, (11) 
i=i 

with ^tr {a Aid Aj) — Sij and the same orthogonality rela- 
tion for party B. 

The optimality of this biorthogonal decomposition 
with respect to the detector efficiencies is shown as fol- 
lows: Consider the most general decomposition W3 = 
bju' A j <8> a Bji where ^'a/Bj are arbitrary (not neces- 
sarily orthogonal) rotated Pauli operators and without 
loss of generality bj > (we can include a minus sign in 
one of the operators). Here, the number of terms M is 
finite and an operator may appear more than once. 

We can express this decomposition in terms of our or- 
thogonal basis o A/Bii 



W 3 



M 

E 

3 



b J a Aj 



•obu (12) 



In Fig. [T]a contour plot of this function is shown: Given 
a certain measured expectation value (W) m , the corre- 
sponding efficiencies rj± ensure that the state is indeed 
entangled. This plot assumes cq + J2 a \ c a\ = lj and can 
be easily redrawn for other decompositions. 

^From eq. © it is obvious that an optimal decom- 
position of such a witness with respect to the needed 
efficiencies rj± is achieved by minimizing ^ Q \c a \. For 
the case of two-qubit witnesses a constructive optimiza- 
tion can be achieved, when arbitrary local Stern-Gcrlach 
measurements (described by Pauli operators or rotations 
thereof) and the identity are allowed. We start with a 
two-qubit witness in its Pauli operator decomposition, 
i.e. 



W = c 00 l 1 + 1® ^2 c 0j' ct J 
/ 3 \ 3 

E cio<Ti ) ® 1 + E c tf <T< ® a i-> ( 10 ) 

\*=1 / i,j=l 

where cri/2/3 = °~xfy/z and c ij S E» The normalization 
condition tr(M^) = 1 leads to coo = \- The three remain- 
ing terms in eq. (|10p can be optimized separately (note 
the special role of the identity) . 

Let us first consider the term W3 = X)i j=i c ij a i ® a j- 
This expression is optimized by doing a singular value 
decomposition of the coefficient matrix (C)y = c^ , i.e. 
UCVi = S, where S is the diagonal matrix that contains 
the singular values s$ . The matrices U and V are orthog- 
onal and have entries Uy and Vij . The new orthogonal ba- 
sis is simply constructed by using the orthonormal rows 



w ith E fe ("fe } ) = E fe (l3k } ) = 1 for all j. The right 
hand sides of eq. (JTTJ) and cq. (fl~2|) are equal. We multiply 
these two expressions by ^ m a Am ® v B m and take the 
trace on both sides. This leads to 



£>tr([ 



CAi ® 0-Bi\ \0~Ar. 



> CfBm}) 



E h a k ] ^^{[dAk ® °Bl] \o Am ® VBm]) ■ (13) 

j,k,l,m 

Orthogonality of the basis is used to get 

= 5>a&>/3&> (14) 

i j,m j 

where we used the fact that the scalar product between 
two normalized vectors (a w and (3 v> ) is less or equal to 
one. This proves the optimality of the decomposition of 
14^3 given in eq. (fTTj) . 

The third term in eq. fit)]) - and, analogously, the sec- 
ond one - can be written as 



W 2 = C10 I E ai ° l 



ciocio 



(15) 



of U and V, i.e. ua% = J2j u ij a j and <*Bi = J2j 



Vij<T n 



where ctio is a rotated Pauli operator, ^ i oq = 1 and 
Sio > 0. Following similar arguments as before, it is easy 
to verify that C10 is optimal. 

The situation for higher dimensions is different: It is 
unfortunately not straightforward to generalize the above 
optimization to higher-dimensional witnesses, because we 
extensively used the fact that a linear combination of 
Pauli operators is again a scaled rotated Pauli operator. 
Also, for multi-partite witnesses the Schmidt decompo- 
sition eq. (|lip does not always exist, such that our opti- 
mization method is not applicable for these cases. 
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FIG. 2: The detector efficiency loophole is closed (assuming 
77+ = 1) if the detector efficiency (77-) is in the hatched area 
for a given measured expectation value {W) m . Here, a de- 
composition with co + ~^2 a \ca\ = 1 is assumed. 



In many experimental situations, e.g. when optical de- 
tectors are used, only the "lost event efficiency" is an 
important issue and the "additional event efficiency" is 
approximately 77+ ~ 1 (l2l . [l3j . This situation further 
simplifies eq. ([8]), and the minimal detector efficiency that 
allows to close the loophole has a simple relation with the 
measured expectation value and the decomposition of the 
witness, namely 



n- > 1 



C + Ea M 



(16) 



This function is shown in fig. where the hatched area 
corresponds to values for which the inequality is fulfilled, 
again assuming cq + ^ a \c a \ = 1. For two qubits the de- 
tection loophole for witnesses can already be closed with 
a detection efficiency 7j_ > |. This bound is sufficient 
e.g. for the optimal two-qubit witness of a Bell state and 
a measured expectation value of (W) m = — \. We con- 



jecture that for r\- 
any witness. 



< 



the loophole cannot be closed for 



III. EXPERIMENTAL EXAMPLES AND 
SUMMARY 



We now want to discuss some prominent experimental 
examples in this context. In ion trap experiments the 
detection efficiencies are close to one (~ 99.8%) [3 

and 

the detection loophole is usually not an issue. Excep- 
tions are many-party-witnesses with expectation values 
close to zero, like the genuine 8-qubit multipartite entan- 
glement witness experiment of Haffner et al. [151 ] . 

Single photon experiments on the other hand are more 
problematic. Using eq. (fit))) we give some explicit exper- 



imental examples for the needed detection efficiencies to 
close the witness loophole: M. Barbieri et al. [l2j im- 
plemented the optimal two-qubit entanglement witness 
to detect a Bell state, where they achieved an expec- 
tation value of —0.493 ± 0.008. In this case the detec- 
tion efficiency needs to be r\- > 0.67. In recent experi- 
ments also multipartite entanglement witnesses were im- 
plemented [131 ]. In this work the three-qubit GHZ en- 
tanglement witness is loophole-free with a detection ef- 
ficiency of 7?_ > 0.91, and the four-partite case needs 
?7_ > 0.94. Single photon detector efficiencies for wave- 
lengths of 700-800 nm are typically around 70% [II], 
such that the global detection efficiency for two qubits 
is circa 50 % , and even lower for more than two subsys- 
tems. The detection efficiencies for multipartite witness 
experiments with photons are thus considerably below 
the needed thresholds. This is a similar situation as for 
the detection loophole in Bell inequalities [4|. However, 
there is a good chance for loophole-free witness experi- 
ments with two qubits, when slightly more efficient de- 
tectors are available. 

In summary, we discussed the detection loophole prob- 
lem for experiments measuring witness operators. As- 
suming the worst case (that may occur due to unknown 
global properties of the detectors) we derived certain in- 
equalities to close such loopholes. These inequalities are 
generally valid for any type of witness operator and de- 
pend on the measured expectation value of the witness, 
its local operator decomposition and the detector effi- 
ciencies. ^From there, detector efficiency thresholds to 
close the loophole are easily calculated. The local decom- 
position of the witness can be optimized such that the 
needed detection efficiencies are minimized. We explic- 
itly presented a constructive optimization for two-qubit 
witnesses. For multi-qubit witnesses the optimal decom- 
position is achieved by minimizing the sum of the ab- 
solute values of the expansion coefficients. In the case 
of higher-dimensional witnesses the optimization is not 
straightforward any more, because it then also depends 
on the type of operator basis. Let us mention that an 
analogous study can be performed, if the witness is de- 
composed into local projectors [l7|; this will be published 
elsewhere. For qubit witnesses we further considered the 
common experimental situation, where additional counts 
can be neglected. Current witness experiments with po- 
larized photons do not close the detection loophole, be- 
cause of the low single photon detector efficiencies. - Fur- 
ther research directions and open problems include the 
optimal local decomposition for higher-dimensional wit- 
nesses, and the case of erroneous detector orientations. 
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